Ancient DNA Research Articles

Archaeozoology Supports a Holistic View on Fish Assessments in Large Rivers—A Case Study from the Volga River: From Quantitative Data and Ancient DNA to Biodiversity Analysis

Knowledge of the historical and present dynamics of populations of migratory fish can promote our understanding of factors affecting their recruitment and abundance. Taxonomic identification of 23,802 bone remains and 13,539 scales of fish from 30 archaeological sites along Volga River revealed that they belonged to 41 different fish species. These data allow for retrospective comparisons and highlight the potential of archaeozoology in conservation biology. Sturgeons and salmonids are vulnerable to the impacts of fishery and climatic change. The sharp decline in the numbers of Starry sturgeon (Acipenser stellatus), Caspian trout (Salmo caspius), and Caspian Inconnu (Stenodus leucichthys) from the Volga in the 17th and 18th centuries was likely related to a cooling period (“Little Ice Age”). At present, the population numbers of all anadromous sturgeons and salmonids of the Volga River are critically low. In the Volga basin over the past two millennia, the sterlet (Acipenser ruthenus) has had a very large population number, high genetic (haplogroups and haplotypes) diversity, and large body sizes. Genetic analysis (aDNA and eDNA) has great potential to expand the knowledge of fish populations along large rivers and to improve long-term biomonitoring. Therefore, analyses of historical data, conventional surveys, as well as the inclusion of genetic approaches complement each other in the development of effective conservation strategies.

Y-STR analysis of highly degraded DNA from skeletal remains over 70years old.

The goal of the following study is to clarify whether the skeletal remains over 70years old from missing persons and their alleged relatives shared identical Y-STR loci. Nowadays, advances in ancient DNA extraction techniques and approaches of using multiple different Y-STRs have significantly increased the possibility of obtaining DNA profiles from highly degraded skeletal remains. Given the ages and conditions of the skeletal remains, ancient DNA extraction methods can be used to maximize the probability of DNA recovery. Considering that information about distant relatives is more relevant for long-term missing persons and alleged family members are male, Y-STR loci analysis is considered the most appropriate and informative approach for determining paternal lineage relationship. In this study, Y-STR genotypes obtained from these alleged relatives were identical to each other and to the alleles of missing persons' consensus profiles at more than 22 loci examined, whilst not being found in Y-STR population database from Y-Chromosome STR Haplotype Reference Database. Therefore, Missing Person No.7 and Missing Person No.18 have a patrilineal relationship with reference samples from Family1 and Family2, respectively. In addition, the fact that Y-STR haplotypes obtained from skeletal remains of missing persons and reference samples are not found in the Han Chinese people from East Asian demonstrates its rarity and further supports a paternal lineage relationship amongst them.

Vegetation and glacier dynamics are sensitive to summer (not winter) warming and the evidence for larch refugia in the ‘Northern Pole of Cold’ inferred from sedimentary ancient DNA and geochemistry

Climate seasonality critically influences the functioning and dynamics of ecosystems in continental areas. The ecological importance of winter temperatures on high-latitude vegetation changes has recently been argued to be largely overlooked in comparison to summer temperatures. The Oymyakon region from eastern Siberia, with its strong continentality of extremely cold winters and moderately warm summers, is ideally suited to study the response of past vegetation to seasonal temperature changes based on long ecological time-series. However, few paleorecords are available from this area. The history of regional glacial activity and potential plant refugia since Marine Isotope Stage (MIS) 3 is not well understood. Here, we present geochemical and plant DNA metabarcoding records retrieved from a sediment core from Lake Ulu in the Oymyakon region, which provides detailed information on glacier and vegetation dynamics over the last 43 cal. ka BP. Our results suggest that glacial fluctuations were primarily driven by summer insolation, and Lake Ulu was likely initiated by glacial retreat during MIS 3. The catchment experienced multiple glacial advance/retreat cycles until the Last Glacial Maximum, and the glaciers fully retreated by 20 cal. ka BP. In addition, a tundra-steppe landscape dominated by Dryas, Papaver, Saliceae, and Anthemideae occupied the catchment for most of the time and began to collapse around 19 cal. ka BP following the expansion of trees and shrubs such as Larix, Betula, Alnus, and Vaccinium. Postglacial plant assemblages in the Oymyakon region exhibit a high sensitivity to summer temperature variations, with minimal impact from winter temperatures. This can be explained by the dominance of summer insolation amplitude, extreme continentality, extended plant growing season, and plant genetic adaptation to cold. Notably, our ancient DNA record show the earliest postglacial expansion of larch in eastern Siberia (around 18.6 cal. ka BP), which is likely related to the presence of local refugia. This implies that the Oymyakon region may be one of the earliest sources for larch recolonization and that more research should be implemented to provide insights into larch expansion and migration, and to better predict the future scenarios for Siberian larch forests.

Ancient DNA indicates a century of overhunting did not reduce genetic diversity in Pacific Walruses (Odobenus rosmarus divergens)

Pacific Walruses (Odobenus rosmarus divergens [Illiger 1815]) are gregarious marine mammals considered to be sentinels of the Arctic because of their dependence on sea ice for feeding, molting, and parturition. Like many other marine mammal species, their population sizes were decimated by historical overhunting in the nineteenth and twentieth centuries. Although they have since been protected from nearly all commercial hunting pressure, they now face rapidly accelerating habitat loss as global warming reduces the extent of summer sea ice in the Arctic. To investigate how genetic variation was impacted by overhunting, we obtained mitochondrial DNA sequences from historic Pacific Walrus samples in Alaska that predate the period of overhunting, as well as from extant populations. We found that genetic variation was unchanged over this period, suggesting Pacific Walruses are resilient to genetic attrition in response to reduced population size, and that this may be related to their high vagility and lack of population structure. Although Pacific Walruses will almost certainly continue to decline in number as the planet warms and summer sea ice is further reduced, they may be less susceptible to the ratcheting effects of inbreeding that typically accompany shrinking populations.

Origin and dispersal history of Hepatitis B virus in Eastern Eurasia

Hepatitis B virus is a globally distributed pathogen and the history of HBV infection in humans predates 10000 years. However, long-term evolutionary history of HBV in Eastern Eurasia remains elusive. We present 34 ancient HBV genomes dating between approximately 5000 to 400 years ago sourced from 17 sites across Eastern Eurasia. Ten sequences have full coverage, and only two sequences have less than 50% coverage. Our results suggest a potential origin of genotypes B and D in Eastern Asia. We observed a higher level of HBV diversity within Eastern Eurasia compared to Western Eurasia between 5000 and 3000 years ago, characterized by the presence of five different genotypes (A, B, C, D, WENBA), underscoring the significance of human migrations and interactions in the spread of HBV. Our results suggest the possibility of a transition from non-recombinant subgenotypes (B1, B5) to recombinant subgenotypes (B2 - B4). This suggests a shift in epidemiological dynamics within Eastern Eurasia over time. Here, our study elucidates the regional origins of prevalent genotypes and shifts in viral subgenotypes over centuries.

Subgenome-aware analyses reveal the genomic consequences of ancient allopolyploid hybridizations throughout the cotton family

Malvaceae comprise some 4,225 species in 243 genera and nine subfamilies and include economically important species, such as cacao, cotton, durian, and jute, with cotton an important model system for studying the domestication of polyploids. Here, we use chromosome-level genome assemblies from representatives of five or six subfamilies (depending on the placement of Ochroma) to differentiate coexisting subgenomes and their evolution during the family's deep history. The results reveal that the allohexaploid Helicteroideae partially derive from an allotetraploid Sterculioideae and also form a component of the allodecaploid Bombacoideae and Malvoideae. The ancestral Malvaceae karyotype consists of 11 protochromosomes. Four subfamilies share a unique reciprocal chromosome translocation, and two other subfamilies share a chromosome fusion. DNA alignments of single-copy nuclear genes do not yield the same relationships as inferred from chromosome structural traits, probably because of genes originating from different ancestral subgenomes. These results illustrate how chromosome-structural data can unravel the evolutionary history of groups with ancient hybrid genomes.

Holme I (Seahenge) and Holme II: ritual responses to climate change in Early Bronze Age Britain

Holme I and II were contemporary, adjacent Early Bronze Age (EBA) oak-timber enclosures exposed intertidally at Holme-next-the-sea, Norfolk, England, in 1998. Holme I enclosed a central upturned tree-stump, its function and intent unknown. Holme II is thought a mortuary structure. Both are proposed here best explained as independent ritual responses to reverse a period of severe climate deterioration recorded before 2049 BC when their timbers were felled. Holme I is thought erected on the summer-solstice, when the cuckoo traditionally stopped singing, departing to the ‘Otherworld’. It replicated the cuckoo’s supposed overwintering quarters: a tree-hole or the ‘bowers of the Otherworld’ represented by the tree-stump, remembered in folklore as ‘penning-the-cuckoo’ where a cuckoo is confined to keep singing and maintain summer. The cuckoo symbolised male-fertility being associated with several Indo-European goddesses of fertility that deified Venus - one previously identified in EBA Britain. Some mortal consorts of these goddesses appear to have been ritually sacrificed at Samhain. Holme II may be an enclosure for the body of one such ‘sacral king’. These hypotheses are considered, using abductive reasoning, as ‘inferences to the best explanations’ from the available evidence. They are supported with environmental data, astronomic and biological evidence, regional folklore, toponymy, and an ethnographic analogy with indigenous Late Iron Age practices that indirect evidence indicates were undertaken in EBA Britain. Cultural and religious continuity is supported by textual sources, the material record and ancient DNA (aDNA) studies.

Patagonian partnerships: the extinct Dusicyon avus and its interaction with prehistoric human communities.

The southern Mendoza province, located in the northern region of Patagonia, was inhabited by hunter-gatherer groups until historic times. Previous archaeological studies have reported canid remains among faunal assemblages, which were assumed to be part of the human diet. However, the taxonomic identification and significance of these canids within human groups have raised questions. In this study, we used ancient DNA analysis, morphological examination and stable isotope analysis (δ13Ccol and δ15N) to re-evaluate the taxonomic assignment of a canid discovered at the Late Holocene burial site of Cañada Seca. Previous morphological identifications suggested that it belonged to the genus Lycalopex, but our results conclusively demonstrate that the individual belongs to the extinct fox species Dusicyon avus. This finding expands Dusicyon avus' known geographical distribution to Patagonia's northern extremity. Furthermore, statistical predictions based on genetic divergence undermine the hypothesis that hybridization between Canis and Dusicyon, facilitated by the introduction of domestic dogs, played a role in the extinction of Dusicyon species. On the other hand, our findings indicate that a Dusicyon avus individual shared a similar diet and was probably buried alongside humans, suggesting a close relationship between the two species during their lives and deaths.

Ancient DNA and osteological analyses of a unique paleo-archive reveal Early Holocene faunal expansion into the Scandinavian Arctic.

Paleo-archives are essential for our understanding of species responses to climate warming, yet such archives are extremely rare in the Arctic. Here, we combine morphological analyses and bulk-bone metabarcoding to investigate a unique chronology of bone deposits sealed in the high-latitude Storsteinhola cave system (68°50' N 16°22' E) in Norway. This deposit dates to a period of climate warming from the end of the Late Glacial [~13 thousand calibrated years before the present (ka cal B.P.)] to the Holocene thermal maximum (~5.6 ka cal B.P.). Paleogenetic analyses allow us to exploit the 1000s of morphologically unidentifiable bone fragments resulting in a high-resolution sequence with 40 different taxa, including species not previously found here. Our record reveals borealization in both the marine and terrestrial environments above the Arctic Circle as a naturally recurring phenomenon in past periods of warming, providing fundamental insights into the ecosystem-wide responses that are ongoing today.

Double feature: First genetic evidence of a mother-daughter double burial in Roman period Austria

A human double burial combined with a horse interment, discovered in 2004 in the eastern cemetery of the Roman town of Ovilava (today Wels/AUT), has been the subject of a comprehensive interdisciplinary study. This burial was remarkable for two reasons: firstly, the arrangement of the two individuals, with individual 1 having an arm around individual 2, and secondly, the presence of equine skeletal remains in the same burial pit directly below the human remains. Based on this finding, an initial assessment suggested a mixed-sex pair buried together with a horse. The horse burial was interpreted as an indicator of an early medieval date. The position of the two human skeletons suggested that a male and female individual had been buried together. A thorough examination of both the human and horse skeletons disproved these initial interpretations. Radiocarbon dating of the human and horse skeletons dated them to the Roman period. In addition, ancient DNA (aDNA) analysis of the human remains corrected the earlier sex assumption, revealing that the individuals were two biological females who were first-degree relatives. The age difference of 15 to 25 years between the two suggests a probable mother-daughter relationship. Thus, the application of scientific methods confirmed a rare, combined human-horse burial from Roman antiquity and established the first genetically documented mother-daughter burial from this period in present-day Austria.

Reconstructing human history from ancient genomes: an interview with Nobel laureate Svante Pääbo

NSR's interview with 2022 Nobel laureate Svante Pääbo to learn about his interesting works in ancient hominine genomes, his perspectives on paleo genome studies and his advice to young researchers.

Comparison and optimization of protocols and whole-genome capture conditions for ancient DNA samples.

Ancient DNA (aDNA)obtained from human remains is typically fragmented and present in relatively low amounts. Here we investigate a set of optimal methods for producing aDNA data by comparing silica-based DNA extraction and aDNA library preparation protocols. We also test the efficiency of whole-genome enrichment (WGC) on ancient human samples by modifying a number of parameter combinations. We find that the Dabney extraction protocol performs significantly better than alternatives. We further observed a positive trend with the BEST library protocol indicating lower clonality. Notably, our results suggest that WGC is effective at retrieving endogenous DNA, particularly from poorly-preserved human samples, by increasing human endogenous proportions by 5x. Thus, aDNA studies will be most likely to benefit from our results.

Reconstructing the Genetic Relationship between Ancient and Present-Day Siberian Populations.

Human populations across a vast area in northern Eurasia, from Fennoscandia to Chukotka, share a distinct genetic component often referred to as the Siberian ancestry. Most enriched in present-day Samoyedic-speaking populations such as Nganasans, its origins and history still remain elusive despite the growing list of ancient and present-day genomes from Siberia. Here, we reanalyze published ancient and present-day Siberian genomes focusing on the Baikal and Yakutia, resolving key questions regarding their genetic history. First, we show a long-term presence of a unique genetic profile in southern Siberia, up to 6,000 yr ago, which distinctly shares a deep ancestral connection with Native Americans. Second, we provide plausible historical models tracing genetic changes in West Baikal and Yakutia in fine resolution. Third, the Middle Neolithic individual from Yakutia, belonging to the Belkachi culture, serves as the best source so far available for the spread of the Siberian ancestry into Fennoscandia and Greenland. These findings shed light on the genetic legacy of the Siberian ancestry and provide insights into the complex interplay between different populations in northern Eurasia throughout history.

Isotopic biographies reveal horse rearing and trading networks in medieval London.

This paper reports a high-resolution isotopic study of medieval horse mobility, revealing their origins and in-life mobility both regionally and internationally. The animals were found in an unusual horse cemetery site found within the City of Westminster, London, England. Enamel strontium, oxygen, and carbon isotope analysis of 15 individuals provides information about likely place of birth, diet, and mobility during the first approximately 5 years of life. Results show that at least seven horses originated outside of Britain in relatively cold climates, potentially in Scandinavia or the Western Alps. Ancient DNA sexing data indicate no consistent sex-specific mobility patterning, although three of the five females came from exceptionally highly radiogenic regions. Another female with low mobility is suggested to be a sedentary broodmare. Our results provide direct and unprecedented evidence for a variety of horse movement and trading practices in the Middle Ages and highlight the importance of international trade in securing high-quality horses for medieval London elites.

X-chromosomal STRs in aDNA kinship analysis.

The analysis of ancient DNA (aDNA) from human skeletal remains can provide useful insights when investigating archaeological finds. One popular application of aDNA is to examine genealogical relationships between individuals recovered at the same archaeological site. For the reconstruction of genealogical relationships, several genetic markers are commonly used: autosomal STRs, mitochondrial lineages (based on SNP-analysis) and Y-chromosomal haplotypes (based on Y-STR-analysis). In this paper, we present the additional opportunities that X-STRs provide in aDNA kinship reconstruction, especially in deficiency cases and for the examination of father-daughter relationships. Possible applications are demonstrated on a range of different kinship reconstructions: confirmation of half-siblingship in the Lichtenstein cave (Germany), exclusion of two potential father-daughter relationships in Goslar (Germany), investigation of three siblingships in Boilstädt (Germany) as well as the confirmation of a father-daughter relationship in Stolpe (Germany). This study shows that the analysis of X-STRs can contribute to the investigation of relationship constellations otherwise difficult to approach (e.g. father-daughter relationships) and that X-STRs are useful to support and complement autosomal STRs, mtDNA and Y-STR data.

Ancient Mitogenomes Reveal the Maternal Genetic History of East Asian Dogs.

Recent studies have suggested that dogs were domesticated during the Last Glacial Maximum (LGM) in Siberia, which contrasts with previous proposed domestication centers (e.g. Europe, the Middle East, and East Asia). Ancient DNA provides a powerful resource for the study of mammalian evolution and has been widely used to understand the genetic history of domestic animals. To understand the maternal genetic history of East Asian dogs, we have made a complete mitogenome dataset of 120 East Asian canids from 38 archaeological sites, including 102 newly sequenced from 12.9 to 1 ka BP (1,000 years before present). The majority (112/119, 94.12%) belonged to haplogroup A, and half of these (55/112, 49.11%) belonged to sub-haplogroup A1b. Most existing mitochondrial haplogroups were present in ancient East Asian dogs. However, mitochondrial lineages in ancient northern dogs (northeastern Eurasia and northern East Asia) were deeper and older than those in southern East Asian dogs. Results suggests that East Asian dogs originated from northeastern Eurasian populations after the LGM, dispersing in two possible directions after domestication. Western Eurasian (Europe and the Middle East) dog maternal ancestries genetically influenced East Asian dogs from approximately 4 ka BP, dramatically increasing after 3 ka BP, and afterwards largely replaced most primary maternal lineages in northern East Asia. Additionally, at least three major mitogenome sub-haplogroups of haplogroup A (A1a, A1b, and A3) reveal at least two major dispersal waves onto the Qinghai-Tibet Plateau in ancient times, indicating eastern (A1b and A3) and western (A1a) Eurasian origins.

Simulation-based benchmarking of ancient haplotype inference for detecting population structure

Paleogenomic data has informed us about the movements, growth, and relationships of ancient populations. It has also given us context for medically relevant adaptations that appear in present-day humans due to introgression from other hominids, and it continues to help us characterize the evolutionary history of humans. However, ancient DNA (aDNA) presents several practical challenges as various factors such as deamination, high fragmentation, environmental contamination of aDNA, and low amounts of recoverable endogenous DNA, make aDNA recovery and analysis more difficult than modern DNA. Most studies with aDNA leverage only SNP data, and only a few studies have made inferences on human demographic history based on haplotype data, possibly because haplotype estimation (or phasing) has not yet been systematically evaluated in the context of aDNA. Here, we evaluate how the unique challenges of aDNA can impact phasing and imputation quality, we also present an aDNA simulation pipeline that integrates multiple existing tools, allowing users to specify features of simulated aDNA and the evolutionary history of the simulated populations. We measured phasing error as a function of aDNA quality and demographic history, and found that low phasing error is achievable even for very ancient individuals (∼ 400 generations in the past) as long as contamination and average coverage are adequate. Our results show that population splits or bottleneck events occurring between the reference and phased populations affect phasing quality, with bottlenecks resultingin the highest average error rates. Finally, we found that using estimated haplotypes, even if not completely accurate, is superior to using the simulated genotype data when reconstructing changes in population structure after population splits between present-day and ancient populations. We also find that the imputation of ancient data before phasing can lead to better phasing quality, even in cases where the reference individuals used for imputation are not representative of the ancient individuals.

Simulation-based benchmarking of ancient haplotype inference for detecting population structure

Paleogenomic data has informed us about the movements, growth, and relationships of ancient populations. It has also given us context for medically relevant adaptations that appear in present-day humans due to introgression from other hominids, and it continues to help us characterize the evolutionary history of humans. However, ancient DNA (aDNA) presents several practical challenges as various factors such as deamination, high fragmentation, environmental contamination of aDNA, and low amounts of recoverable endogenous DNA, make aDNA recovery and analysis more difficult than modern DNA. Most studies with aDNA leverage only SNP data, and only a few studies have made inferences on human demographic history based on haplotype data, possibly because haplotype estimation (or phasing) has not yet been systematically evaluated in the context of aDNA. Here, we evaluate how the unique challenges of aDNA can impact phasing and imputation quality, we also present an aDNA simulation pipeline that integrates multiple existing tools, allowing users to specify features of simulated aDNA and the evolutionary history of the simulated populations. We measured phasing error as a function of aDNA quality and demographic history, and found that low phasing error is achievable even for very ancient individuals (∼ 400 generations in the past) as long as contamination and average coverage are adequate. Our results show that population splits or bottleneck events occurring between the reference and phased populations affect phasing quality, with bottlenecks resulting in the highest average error rates. Finally, we found that using estimated haplotypes, even if not completely accurate, is superior to using the simulated genotype data when reconstructing changes in population structure after population splits between present-day and ancient populations. We also find that the imputation of ancient data before phasing can lead to better phasing quality, even in cases where the reference individuals used for imputation are not representative of the ancient individuals.

Exploring the complex pre-adaptations of invasive plants to anthropogenic disturbance: a call for integration of archaeobotanical approaches.

Pre-adaptation to anthropogenic disturbance is broadly considered key for plant invasion success. Nevertheless, empirical evidence remains scarce and fragmentary, given the multifaceted nature of anthropogenic disturbance itself and the complexity of other evolutionary forces shaping the (epi)-genomes of recent native and invasive plant populations. Here, we review and critically revisit the existing theory and empirical evidence in the field of evolutionary ecology and highlight novel integrative research avenues that work at the interface with archaeology to solve open questions. The approaches suggested so far focus on contemporary plant populations, although their genomes have rapidly changed since their initial introduction in response to numerous selective and stochastic forces. We elaborate that a role of pre-adaptation to anthropogenic disturbance in plant invasion success should thus additionally be validated based on the analyses of archaeobotanical remains. Such materials, in the light of detailed knowledge on past human societies could highlight fine-scale differences in the type and timing of past disturbances. We propose a combination of archaeobotanical, ancient DNA and morphometric analyses of plant macro- and microremains to assess past community composition, and species' functional traits to unravel the timing of adaptation processes, their drivers and their long-term consequences for invasive species. Although such methodologies have proven to be feasible for numerous crop plants, they have not been yet applied to wild invasive species, which opens a wide array of insights into their evolution.

Mitochondrial genomes of Pleistocene megafauna retrieved from recent sediment layers of two Siberian lakes

Ancient environmental DNA (aeDNA) from lake sediments has yielded remarkable insights for the reconstruction of past ecosystems, including suggestions of late survival of extinct species. However, translocation and lateral inflow of DNA in sediments can potentially distort the stratigraphic signal of the DNA. Using three different approaches on two short lake sediment cores of the Yamal peninsula, West Siberia, with ages spanning only the past hundreds of years, we detect DNA and identified mitochondrial genomes of multiple mammoth and woolly rhinoceros individuals—both species that have been extinct for thousands of years on the mainland. The occurrence of clearly identifiable aeDNA of extinct Pleistocene megafauna (e.g. >400 K reads in one core) throughout these two short subsurface cores, along with specificities of sedimentology and dating, confirm that processes acting on regional scales, such as extensive permafrost thawing, can influence the aeDNA record and should be accounted for in aeDNA paleoecology.